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Abstract 

We predict and tr ansverse momentum distributions from proton- 
proton and nuclear collisions at RHIC and LHC. A resummation formalism 
with power corrections to the renormalization group equations is used. The 
dependence of the resummed QCD results on the non-perturbative input is 
very weak for the systems considered. Shadowing effects are discussed and 
found to be unimportant at RHIC, but important for LHC. We study the 
enhancement of power corrections due to multiple scattering in nuclear col- 
lisions and numerically illustrate the weak effects of the dependence on the 
nuclear mass. 
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I. INTRODUCTION 



With the commissioning of the Relativistic Heavy Ion CoUider (RHIC), nuclear coUision 
physics entered the coUider era. New phenomena have been presented in Au + Au colhsions 
at the center-of-mass energy y/s = 130 GeV/nucleon Negative hadron multiphcities 
are significantly increased compared to proton-antiproton collisions 0. Some of the data, 
such as the antiproton-to-proton ratios ||^ , are reported to contradict perturbative Quantum 
Chromodynamics (pQCD) calculations. Initial efforts to understand hadron spectra from 
nuclear reactions at RHIC using pQCD have been published recently However, there 

is still a strong debate on how to apply pQCD to hadron production even in proton-proton 
(pp) collisions. With RHIC running at the full energy, and experiments at the Large Hadron 
Collider (LHC) at several TeV/nucleon soon to follow, "clean" processes, where pQCD works 
well in pp collisions, are needed to test pQCD in nuclear reactions at coUider energies. It is 
particularly important to distinguish between expected nuclear effects and "new physics" . 

The production of and Z° bosons has been extensively studied at Tevatron energy 
{y/s = 1.8 TeV). Perturbative QCD proved to be very successful in explaining the CDF 
and DO data. Due to the large mass of and Z^, and no final state rescattering in 
their production process, power corrections are expected to be small. Therefore, and 
Z° production could provide a bench mark test for pQCD at RHIC and LHC in both pp and 
nuclear collisions. Luminosity remains a concern at RHIC, in particular in pp collisions. 

The bulk of the data for and production is concentrated in the small transverse 
momentum (pt) region, where Pt is much smaller than the corresponding heavy boson mass 
(M). When pj- -C M, the pr distributions calculated order- by-order in in conventional 
fixed-order perturbation theory receive a large logarithm, ln(M^/p^), at every power of as, 
even in the leading order in as- Therefore, at sufficiently small pt, the convergence of the 
conventional perturbative expansion in powers of as is impaired, and the logarithms must 
be resummed. 

Resummation of the large logarithms in QCD can be carried out either in p^-space 
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directly, or in the so-called "impact parameter", fe-space, which is a Fourier conjugate of 
the pT-space. (Since we are also interested in nuclear collisions, we reserve the notation 
b for the usual impact parameter associated with the geometry of colliding heavy nuclei, 
and denote the mathematical tool used in the resummation procedure by the symbol b.) 
Using the renormalization group equation technique, Collins and Soper improved the 6-space 
resummation to resum all logarithms as singular as, ln™'(M^/p^)/p|., as ^ Q]. In the 
frame work of this renormalization group improved 6-space resummation, Collins, Soper, 
and Sterman (CSS) derived a formalism for the transverse momentum distribution of vector 
boson production in hadronic collisions In the CSS formalism, non-perturbative input is 
needed for the large b region. The dependence of the pQCD results on the non-perturbative 
input is not weak if the original extrapolation proposed by CSS is used. Recently, a new 
extrapolation scheme was proposed, based on solving the renormalization group equations, 
including power corrections [jlO[- Using the new extrapolation formula, the dependence of 
the pQCD result on the non-perturbative input was significantly reduced. This was achieved 
without matching the fixed-order calculations. The results agree with the DO nd CDF data 
very well in the entire pt interval from px ^ 1 GeV to px as large as the the vector 
mass. In addition, in this new extrapolation formalism, power corrections can be studied 
systematically for pp to proton-nucleus (pA) to nucleus-nucleus (AB) collisions. 

Factorization plays a key role in pQCD predictions of observables. At RHIC and LHC, 
nuclear collisions reach an unprecedented energy scale. Factorization can not be extrapo- 



lated easily to AB collisions from hadron-hadron collisions [jTT| , and power corrections may 
become important [T^. There are two types of power corrections: (1) Power corrections 
directly to the physical observables (type-I) and (2) power corrections to the evolution of 
renormalization group equations (type-II). Type-I power correction are proportional to pow- 
ers of (Aqcd/Q), Q being the physical large scale (momentum). These are small for 
and production as a result of the large mass of the particles in question. Type-II power 
corrections are proportional to powers of (Aqcd//^), with evolution scale fi. Therefore, 
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physical observables carry the effect of type-II power corrections for all ^[Qq,Q]^ with the 
boundary condition at the scale Qq. Even with large mass, and can still carry a large 
effect of type-II power corrections. However, Ref. |10[ provides a way to study the effect of 



type-II power corrections numerically. Based on the numerical results presented in Ref. [10 



the effect of type-II power corrections is also expected to be small for heavy vector boson 
production in pp collisions. In nuclear collisions, type-II power corrections will be enhanced 
by the nuclear size. Here we study the effects of enhanced power correction in pA and AB 
collisions. 

The rapidity distribution of W"^ and bosons at LHC was calculated recently and 
the effect of shadowing was discussed in this context |13[. The decay spectrum of Z^ was 
suggested as an alternative reference process for J/ip suppression at LHC, since the Drell- 
Yan lepton pairs at low mass are expected to be dominated by hh decay. As most of the 
and Z^ bosons observed in the experiment will be in the small transverse momentum 
region, it is of great interest to study the nuclear effects on massive gauge boson transverse 
momentum distributions. These effects include both the nuclear modification of the parton 
distribution function (shadowing), and the initial state power corrections. The second part 
of the present paper is devoted to these subjects. 

The rest of the paper is organized as follows. In Section ||, we review the CSS 6-space 
resummation formalism and the new procedure |]10| to extrapolate the pQCD calculation to 
the large h region is outlined. We give predictions for and Z^ transverse momentum 
distributions from pp collision at both RHIC and LHC energies. The sensitivity of the 
prediction to the non-perturbative parameters is discussed. Nuclear effects on and Z° 



transverse momentum distributions in pA and AB collisions are studied in Sections |ITl| and 
P^ . In Section [V|, we present our conclusions. To make the paper reasonably self-contained, 
the Appendix summarizes relevant details of the CSS formalism. We use h = c = 1 units. 
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II. TRANSVERSE MOMENTUM DISTRIBUTION OF HEAVY VECTOR 
BOSONS IN PROTON-PROTON COLLISIONS 



In this Section, we briefly review the CSS 6-space resummation formahsm (Section p.1 A| ). 
In Section [il B| we introduce the recently-proposed new formalism to extrapolate the per- 



turbative calculation to the large b region [1^], and discuss the transverse momentum dis- 
tribution of heavy vector bosons in pp collisions at RHIC and LHC energies. 

A. Collins-Soper-Sterman (CSS) formalism 

In the production of heavy vector bosons (of mass M) , when p^ <^ M, the Pt distribution 
calculated order-by-order in in conventional fixed-order perturbation theory receives a 
large logarithm, ln(M^/p^), at every power of as- Even at the leading order in as, the cross 
section da/dM'^ dydp\ contains a term proportional to (os/Pt) ln(M^/py) coming from the 
partonic subprocess q + q —* V + g, where V = •j*, or Z^. Beyond the leading order, 
we actually get two powers of the logarithm for every power of as, due to soft and collinear 
gluons emitted by the incoming partons. Therefore, at sufficiently small pt, the convergence 
of the conventional perturbative expansion in powers of as is impaired, and the logarithms 
must be resummed. 

For vector boson production in hadronic collisions between hadrons Ha and Hb, hA+hs — »• 
V{M) + X with V = 7*, W^, or Z^, the CSS resummation formalism yields the following 
generic form 

where xa = Mj^ and xb = e'^M/i/i, with rapidity y and collision energy ^/s. In 
Eq. (|l|), the W term dominates the px distributions when px <^ M, and the Y term gives 
corrections that are negligible for small pt, but become important when p^ ~ M. In the 
CSS formalism, the function W appears as a superposition. 
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W{b, M,xa,xb)=Y. xa,xb) cx,j^v{M) , (2) 

where aij^v{M) is the lowest order cross section for a quark-antiquark pair of invariant 
mass M to annihilate into a vector boson V, and the X^ij runs over all possible quark and 
antiquark flavors that can annihilate into the vector boson at the Born level 0. In Eq. @, 
Wij{b, M,xa,xb) is an effective flux to have partons of flavor i and j from the respective 
hadrons Ha and hs, and it has the form 

Wij{b, M, XA, xb) = e^(^'*') Wij{b, c/b, XA, xb) , (3) 

where S{b,M) is given in the Appendix, and c is a constant of order unity The 
functions Wij{b,c/b,XA,XB) in Eq. (H) depend on only one momentum scale, 1/6, and 
are perturbatively calculable as long as 1/6 is large enough. All large logarithms from 
ln(l/6^) to In(M^) in Wij{b, M,xa,xb) are completely resummed into the exponential fac- 
tor exp[S{b,M)]. 

Since the perturbatively resummed Wij{b, M,xa,xb) in Eq. (j^) is only reliable for the 
small b region, an extrapolation to the large b region is necessary in order to complete 
the Fourier transform in Eq. ([1|). This is the point where the dependence of the result on 
the nonperturbative parameters enters into the formalism, possibly limiting the predicative 
power of pQCD , especially for small px- We will next discuss this in detail. 



B. Improved extrapolation formula with power corrections 

The CSS formalism in its original form introduces a modification to the output of the 



perturbative calculation (see Appendix, Eq. (p5D). The size of the modification depends 
on the non-perturbative input parameters, and can be as large as 50% of the perturbative 
result. Therefore the predictive power of the CSS formalism was questioned and new efforts 
have been devoted to resum the large logarithms directly in pt space p4| , p3| . 

In oder to improve the predictive power of the CSS 6-space resummation formalism, 
the following new form was proposed [|lOl by solving the renormalization group equation 
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including power corrections: 



W(b, M,xa,xb) 



W{b,M,XA,XB) 



h<hr. 



W{bmax, M, xa, xb) F^^{b, M, Xa, Xb] bmax, «) b > b„ 
where the nonperturbative function in Eq. (H) is given by 



(4) 



F^^{b, M, xa,xb; bmax, a) = exp <{ - In 



M'b: 



2U2 



9i{ibr-{bLxr)+92ib'-b. 



;2 l2 

max 



92 - b'm 



(5) 



In Eq. (^) bmax is a parameter to separate the perturbatively calculated part from the 
non-perturbative input. Unlike in the original CSS formalism, when b < bmax, the perturba- 
tively calculated W{b, M,xa,xb) is not altered and is independent of the nonperturbative 
parameters. In addition, the 6-dependence in Eq. is separated according to different 
physics origins. The (6^)"-dependence mimics the summation of the perturbatively cal- 
culable leading power contributions to the kernels K and G in Eq.s (|T3|) and (0) of the 
Appendix, respectively, to all orders in the running coupling constant q;s(/u) with the scale 
/i running into the nonperturbative region []T^ . The 6^-dependence of the g2 term is a direct 
consequence of dynamical power corrections to the renormalization group equations of the 
kernels K and G, and has an explicit dependence on M. The g2 term represents the effect 
of the non-vanishing intrinsic parton transverse momentum. As we discuss later, these two 
terms behave differently in nuclear collisions. 

A remarkable feature of the 6-space resummation formalism is that the resummed ex- 
ponential factor exp[S'(6, M)] (see Eq. (^) suppresses the 6-integral when b is larger than 
1/M. Therefore, when M ^ fiQ, (where /iq ~ ^/bmax represents the scale at which pQCD 
starts to be applicable), it is possible that the Fourier transform in Eq. ([^) is dominated 
by a region of b much smaller than l//io, and the calculated pt distributions are insensi- 
tive to the nonperturbative information at large b. In fact, it was shown using the saddle 
point method that, for a large enough M, a QCD perturbation calculation is valid even 
at Pt = [p!7||9|. As discussed in Ref. |T^, the value of the saddle point strongly depends 



d%e'P^^-'Wib,M,XA,XB) = — / dbbJoipTb)W{b,M,XA,XB) , (6) 

/TT JO 



on the collision energy ^/s, in addition to its well-known dependence. Because of the 
steep evolution of parton distributions at small x, the y/s dependence of the W in Eq. (|l]) 
significantly decreases the value of the saddle point and improves the predictive power of 
the 6-space resummation formalism at collider energies. 

To display the saddle point clearly, let us rewrite the first term on the right-hand side of 
Eq. ([l|), taking into account that there is no preferred transverse direction and that W in 
Eq. (|l]) is a function of 6 = |6| only: 

1 _ ^ _ _ 1 POO 

— — / d'b e^P^-" W{b, M, XA, ^b) = ^ " ' 

where W contains the exponential factor exp[S'(6, M)], and Jq{z) with z = prb represents 
the Bessel function of order zero. When px > 0, the Bessel function suppresses the large 
b region of the integration. Because the argument of the Bessel function is proportional to 
Pt, the large b region is more suppressed if pt is larger. In the following, we focus on the 
saddle point aX pr = 0. Fig. 1 is a graph (on an arbitrary scale) of the integrand of the 
right hand side of Eq. (|^) with pr = for the case of production at the collision energies 
where we intend to carry out our calculations: y/s = 350 GeV (dotted line), i/i = 500 GeV 
(short dashes), and at = 5.5 TeV (solid). The position of the saddle point decreases 
as the collision energy increases. We will show that the predictive power of the formalism 
is very good at RHIC. As a consequence of the behavior of the saddle point position, we 
expect excellent precision at y/s = 5.5 TeV. At the LHC pp energy of ^/s = 14 TeV, the 
saddle point moves to even smaller values of b, improving the accuracy of the calculation 
even further. 

The parameters gi and a of Eq. (^) are fixed by the requirement of continuity of the 
function W{b) and its derivative at 6 = bmax- (The results are insensitive to changes of bmax 



in the interval 0.3 GeV-^ < bmax ^ 0.7 GeV'^ We use bmax = 0.5 GeV-^ in Ref. and 
here.) The value of g2 and g2 can be obtained by fitting the low-energy Drell-Yan data, and 
are taken to be g2 ~ 0.27 GeV^ and g2 ~ 0.01 GeV^ following Ref.s 110] and |TB|. As the b 



dependence of the g2 and g2 terms in Eq. (^ is identical, it is convenient to combine these 
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terms and define 



G2 = ln( ^)<72 + 92 . (7) 



C2 



Using the values of tfie parameters listed above, we get G2 ~ 0.4 GeV^ for and pro- 
duction in pp collisions. The parameter G2 can be considered the only free parameter in the 
non-perturbative input in Eq. (^), arising from the power corrections in the renormalization 
group equation. Our results are not sensitive to small variations of G2 around its estimated 
value. We use G2 = 0.4 GeV^ in the following. An impression about the importance of 
power corrections can be obtained by comparing results with G2 = 0.4 GeV^ to those with 
power corrections turned off {G2 = 0). We therefore define the ratio of these quantities: 

The cross sections in the above equation and in the results presented in this paper have 
been integrated over rapidity and invariant mass squared from Eq. (|l]). For the integration 
over dM'^, we use the "narrow width approximation" |T9|. With respect to rapidity, we 
integrate over the interval [—1, 1] for RHIC, and over [-2.4,2.4] for LHC, representing the 
central acceptances of the appropriate detectors. For the parton distribution functions, we 



use the CTEQ5M set ||2^ in the present work. 

The following figures display differential cross sections and the corresponding ratios 
for and production as functions of px, in order of increasing energy. (In what follows, 
by production cross section we mean the cross section for + W~.) In Fig. 2(a) we 
show the differential cross section for Z^ production in pp collisions at a/s = 500 GeV, 
where a RHIC pp run is expected. The accompanying Fig. 2(b) displays of Eq. (H). It 
can be seen that the ratio converges to unity (i.e. very small power corrections) with 
increasing px, and the effect of power corrections is about five percent for p^ = 0. However, 
when pt > 2 GeV, the deviation of from unity is under two percent. A very similar 
conclusion can be drawn from Fig. 3, concerning production in pp collisions at the same 
energy: the dependence on the nonperturbative input is weak even at the lowest transverse 
momenta. 
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Fig.-s 4 and 5 contain similar information for and production in pp collisions at 
the planned LHC nucleus-nucleus collision energy of ^/s = 5.5 TeV, while the cross section 
and are given for the LHC pp energy of y/s = 14 TeV in Fig.-s 6 and 7. At all energies, 
the magnitude of the production cross section is about 3 to 5 times that of the Z° 
production cross section. Otherwise, the results for and production are very similar, 
because the saddle points are almost at the same place in b space in the two cases. As it has 
been pointed out in Ref flOl, the position of the saddle point is determined by two terms: 
the first term is inversely proportional to the mass of the vector boson, and the second term 
is proportional to the derivative of the parton distribution function with respect to the scale 
at xa or xb- At the collider energies considered here, the second term is negative. While 
M]Y is a little smaller than Mz, and thus the first term is larger for than for Z^, the 
contribution from the second term will effect the saddle point for more than that of the 
Z^ with a small xa and xb- 

As expected from Fig. 1, the dependence on the nonperturbative input decreases with 
increasing collision energy. The ratio is smaller than one percent at LHC for both 
^/s = 5.5 TeV and ^/s = 14 TeV, even when pt = 0. These results confirm that the 
predictive power of the formalism is very good at RHIC and excellent at LHC. They also 
show that the effect of power corrections is very small at LHC for almost the whole Pt region 
and it is also negligible for RHIC when px is larger than 2 GeV. 



III. ISOSPIN AND SHADOWING EFFECTS IN NUCLEAR COLLISIONS 

In lack of nuclear effects on the hard collision, the production of heavy vector bosons in 
nucleus-nucleus (AB) collisions should scale, compared to the production in pp collisions, 
as the number of hard collisions, AB. However, there are several nuclear effects on the 
hard collision in a heavy-ion reaction. We distinguish three categories of these effects: 
isospin effects, the modification of the parton distribution function in the nucleus, and the 
enhancement of power corrections by the nuclear size. In this Section we discuss the first 
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two of these nuclear effects. Power corrections will be treated separately in the next Section. 

As the parton distribution of neutrons is different from that of the protons, the production 
cross section of heavy vector bosons in proton-neutron and neutron-neutron interactions 
differs from the corresponding production cross section in pp collisions. This difference is 
the source of the so-called isospin effects. To study the isospin effects, let us introduce 

da{pT,ZA/A,ZB/B) / da{pT) , . 

'^-''^^^^ " dfr / ' 

where Za and Zb are the atomic numbers and A and B are the mass numbers of the colliding 
nuclei, and the cross section da{pT, Za/A, Zb/ B)/dpT has been averaged over AB, while 
da{pT) / dpT is the pp cross section (with 6*2 = 0). Since the uu coupling in Z^ production 
is replaced by ud coupling to produce bosons, we expect Riso to be complementary for 
and Z^ production, i.e. if one of these ratios is larger than one, the other ratio should 
be smaller than unity. 

Next we turn to the phenomenon of shadowing. We use the term shadowing in the general 
sense, referring to all modifications of the parton distribution function in the environment 
of the nucleus |jT^. This includes, in different regions of x, the phenomena of shadowing (in 
the strict sense), anti-shadowing, and the EMC effect. It is also important to note that we 
define shadowing effects as leading twist effects, not including power corrections |]TT1. Thus, 
shadowing effects and power correction are clearly separated in the present work. 

In general, shadowing is expected to be a function of x, the scale Q, and of the position 
in the nucleus. The latter dependence means that in heavy-ion collisions, shadowing should 
be impact parameter (6) dependent. The parameterizations of shadowing in the literature 
take into account some of these effects, but no complete parameterization exists to date 
to our knowledge. For example, the HIJING parameterization includes impact parameter 
dependence, but does not deal with the scale dependence pT| , p2| . On the other hand, the 



EKS98 ||2^ and HKM ||2^ parameterizations have a scale dependence, but do not consider 
impact parameter dependence. (The latter parameterizations have been compared recently 
p5|.) Since in this paper we concentrate on impact-parameter integrated results, where 
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the effect of the 6-dependence of shadowing is relatively unimportant 0, and we focus 
more attention on scale dependence, we find the EKS98 shadowing parameterization most 



Thus, Rsh incorporates the effects of both shadowing and the isospin composition of the 
nuclei. Furthermore, in this Section we take G2 = throughout, switching off power correc- 
tions. The enhanced power corrections in pA and AB collisions will be discussed separately 
in the next Section. 

Fig. 8 displays Riso (dashed lines) and Rsh (full lines) as functions of pt for production 
(Fig. 8(a)) and production (Fig. 8(b)), respectively, at y/s = 350 MeV. The same 
information is presented at ^/s = 5.5 TeV in Fig. 9. We find that Riso > 1 for production 
and Riso < 1 for production at both energies, complementing each other, as expected. 
How close Riso is to unity depends on the region of x that has the dominant contribution to 
the cross section. At LHC, x ~ 0.02, and the magnitude of the isospin effects is about 2%. 
This is because when x is in this range, the difference of the distributions of u{u) and d{d) 
quarks in the nucleus is very small. At RHIC, where the dominant x ~ 0.26, the isospin 
effects are about 20% in the small pt region. We have also checked that the scale dependence 
of isospin effects is not important. In Fig. 8, we also show a dotted line, which is the ratio 
of Rsh to Riso, and thus can be considered to display the effect of "pure shadowing" . Pure 
shadowing appears as an approximately constant reduction of 3-5% at x ~ 0.26. 

The appearance of Rsh (full lines) is more surprising at y/s = 5.5 TeV (Fig. 9). Since 
at LHC, even at pt = 90 GeV, x ~ 0.05, we are still in the "strict shadowing" region. 
Therefore, the fact that Rsh > 1 for 20 GeV ^ ^ 70 GeV is not "antishadowing" ; rather, 
it is a consequence of the change of the shape of the cross section from pp to AB reactions. 
In pp collisions, the maximum of the cross section is at ~ 3.7 GeV. In AB collisions, the 
peak is shifted by about 0.4 GeV in the direction of large transverse momenta. On the other 



appropriate for our present purposes. Therefore we use EKS98 shadowing in this work. 
Similarly to Eq. (^ for isospin effects, we define a shadowing ratio. 
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hand, the contributions from the scale ~ j9t to M have been resummed. Resummation 
plays a more important role for small px than for large pt- In EKS98, the strict shadowing 
effect is more important for small scales, such as pt ~ 2 — 5 GeV, than for px ~ M. Thus, 
a relatively small shift in the peak position brings about the rise of Rsh above unity. The 
shadowing ratio Rsh dips below one again at a certain pt, to approach the fixed-order pQCD 
result when px ~ M. The situation in Fig. 9 can be understood by noting that shadowing in 
the resummed transverse momentum distributions is sensitive to the scale used to describe 
the nuclear effects on the parton distribution function. 

After this discussion one may want to revisit Fig. 8 and ask why is Rsh flat at RHIC 
energies. The answer lies in the observation that the EMC region of x dominates the cross 
section at i/i = 350 MeV. There is no sensitivity to the scale in that region in EKS98. 

In summary, the shadowing effects in the pt distribution of heavy bosons at RHIC and 
LHC are sensitive to the scale of the nuclear parton distribution. The available data on the 



scale dependence of nuclear parton distribution functions are very limited |^3[. Theoretical 
studies (such as EKS98) are based on the the assumption that the nuclear parton distribution 
functions differ from the parton distributions in the free proton, but obey the same DGLAP 
evolution P3[. Therefore, the tranverse momentum distribution of heavy bosons at LHC 
in Pb + Pb collisions can provide a further test of our understanding of the nuclear parton 
distributions. 



IV. ENHANCED POWER CORRECTIONS IN NUCLEAR COLLISIONS 

We have seen in Section that the effect of power corrections on the resummed cross 
section is very weak in pp collisions at both RHIC and LHC. In pA and AB collisions, power 
corrections will be enhanced due to rescattering in the nucleus. In this Section we discuss 
the enhancement of power corrections in nuclear collisions. 

As g2 represents the effect of the partons' non- vanishing intrinsic transverse momentum, 
it should not have a strong nuclear dependence. On the other hand, g2 arises as a result of 
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dynamical power corrections, and should be enhanced by the nuclear size, i.e. proportional 
to A^^^. Considering the enhancement, we get G2 ~ 0.8 GeV^ for p + Au and G2 ~ 1.5 
GeV^ for Pb + Pb reactions. 

In order to study the effects of power corrections, we now define Rq"^ in nuclear collision 
as follows: 



where the numerator represents the result of the full calculation, with isospin effects, shad- 
owing, and power corrections taken into account, and the denominator was introduced in 



Fig.lO(a) and Fig.ll(a) present the da^'^^''''^\pT, Za/A, Zb/ B) for production at 
a/s = 350 GeV in p + Au reactions and y/s = 5.5 TeV in Pb + Pb collisions, respectively. 
In Fig. 10(b) the solid line is Rq'^ in p + Au collisions and the dashed line is Rg2 in PP 
collisions at the same energy. At ~ GeV, R^^ is about 0.88, in contrast to about 0.95 
in PP collisions. The effects of power corrections are enhanced about three times in p + Au 
collision relative to pp collisions. However, even for ~ 0, the effect of power corrections is 
only about 10%, and for px > 2 GeV, the effect remains under 5%. Thus the effect of power 
corrections is weak in p + Au collisions at RHIC. As for the situation at LHC in Pb + Pb 
collisions, displayed in Fig. 11, the effects of power corrections appear to be enhanced by a 
similar factor (about three) from pp to Pb + Pb collisions at LHC. The reason for the similar 
enhancement is the interplay of the higher energy and the larger value of G2 mentioned 
above. In any case, the enhanced power corrections remain under 3% in the [0,80] GeV pt 
interval for Pb + Pb at LHC. The results for Z^ production are similar to those for 
production. 



R\ 



da^^^''^\pT,ZA 
dpT 




(11) 



connection with Eq. (|T0|) and contains no power corrections {G2 = 0). 



V. DISCUSSION AND CONCLUSION 



We have demonstrated that the predictive power of a resummed pQCD calculation for 
heavy boson production is very good at RHIC and excellent at LHC, because the depen- 
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dence of the cross sections on the non-perturbative input is weak. This makes and 
production rehably calculable at the energies considered. On the experimental side, 
and Z° production should be an early observable at LHC (with sufficiently high production 
rates) in both pp and nuclear collisions. At RHIC, production rates may remain on the 
low side, in particular in pp collisions. In all cases, the cross sections peak at relatively low 
transverse momenta, and we calculate the peak position for all systems considered. 

We have examined the effects of the isospin composition of nuclei and of shadowing 
in nuclear coUisions. We find that isospin effects are modest in and Z° production 
at the energies considered, while shadowing is relatively unimportant at RHIC, but plays 
an important role at LHC in shaping the transverse momentum distributions. This comes 
about because of the sensitivity of shadowing effects to the scale of the parton distribution 
functions. We also studied the magnitude of power corrections in pp and nuclear collisions. 
Power corrections are small in pp collisions, and, even though part of these corrections is 
enhanced in nuclear coUisions by a factor of A^/^, they remain below 10% (even at very low 
transverse momenta) in AB reactions. 

As indicated by the above observations, the resummed pQCD formalism yields accurate 
results for and Z° transverse momentum distributions at RHIC and LHC. The produc- 
tion of heavy bosons can therefore serve as a further test of the nuclear parton distributions. 
In more general terms, and Z° production at RHIC and LHC provides a bench mark 
test of resummed pQCD for pp and nuclear collisions at these energies. We hope that our 
predictions will soon be compared to data. 
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VI. APPENDIX 



Here we summarize relevant details of the CSS formalism concerning the structure of 
the terms in Eq. We address, in particular, the construction of the exponential factor 
exp[S{b,M)] in Eq. (|). 

It was shown in Ref. that for b <^ I/Aqcd, the function W{b, M,xa,xb) of Eq. (|I]) 
is directly related to the singular parts of the pt distribution as pt — * 0. More precisely, 
W{b, M,xa,xb) includes all singular terms like S'^{pt) and [\n^ {p"^ / p"^) / p'^jreg with m > 0. 
Less singular terms as ^ are included in the Y term of Eq. (|T|). The QCD resummation 
of the large logarithms in the CSS formalism is achieved by solving the evolution equation 
for the Wij of Eq. (|), 

d ~ ^ r ~ ... 

-g^^^ijib, M, XA, xb) = [K{bfi, + G(M//i, a,(/i))J W,,ib, M, xa, xb) , (12) 

and the corresponding renormalization group equations for the kernels K and G, 

d ~ 1 

K{bn,asiiJ,)) = --7i^(as(/i)) , (13) 



d 1 

G{M/fi, a,(/i)) = -7x(as(/u)) . (14) 



The anomalous dimension 7/^(as(yu)) = J2n=ilK\'^s{l^) /'n')^ in Eqs. ([TB| ) and (P^ ) is per 



turbatively calculable The renormalization group equations (^) and (H) for K and G 
ensure the correct renormalization scale dependence, 

^»'(M/,..,.b)=0. (15) 

The solution given in Eq. (^ corresponds to simultaneously solving the evolution equation 
(|l^) from ln(c^/62) to ln(M2), and the renormalization group equations (|13D and ([T^ ) from 



ln(c^/62) to ln(/i2) and from ln(M2) to ln(/i2), respectively, where c is a constant of order 
unity. 

Integrating Eq. (|13D from \n{c^/b'^) to \n{ij?), and Eq. (|14[) from ln(M2) to \n{ij^), one 
derives 
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K{bfi, + G'(M//i, 



(16) 



where A is a function of 7x(Q;s(/i)) and K{c,as{jj)) while -B depends on both K{c,as{M)) 
and G{1, as{M)). The functions A and B do not have large logarithms and have perturbative 
expansions A = En=i^^"^(as/7r)" and B = E„=i 5(")(a,/7r)'^, respectively. The first two 
coefficients in the perturbative expansions are known P,|26[|: 



^(1) 


— Cp , 




_ Cf ' 
2 






5(2) 


^(1 



'67_ TT^^ 

18 6 



10 



IT 



7 - 12C(3) 



2 



11 



-TT 



193 



+ 3C(3) 



17 2 2 

TT 

6 9 



(17) 



where = 3 for SU(3) color, the color factor Cp = {N^-1)/2N = 4/3 for quarks, Tr = 1/2, 
and is the number of active quark flavors. The coefficients A and B given in Eq. (|1^ are 
derived from the general expressions in Ref. P] with renormalization constant c = 2e~'''^, 
where •je ~ 0.577 is Euler's constant p. 

Substituting Eq. (0) into Eq. (|13), and integrating from ln(c^/6^) to In(M^), one obtains 



the Wij given in Eq. (|^) with 



S(b, M) 



Af2 



T.I 



c2 y? 



/MM 

In ( — j v4(a,(^)) + B{as{fi)) 



(18) 



In Eq. (^, all large logarithms from ln(c^/6^) to In(M^) in Wij{b, M, xa, xb) are completely 
resummed into the exponential factor exp[S'(6, M)], leaving the Wij{b, c/b, xa,xb) with only 
one momentum scale 1/6. The Wij{b, c/b, xa, xb) in Eq. (^ is then perturbatively calculable 
when the momentum scale 1/6 is large enough, and is given by 



Wij{b, c/b, Xa, Xb) = fi/A{xA, At = c/b) fj/B{xB, yU = c/b) . 
The functions fi/A and fj/B are the modifled parton distributions PJT4 
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(19) 



fi/A{xA,f^) = ^C'»/a(a;A/^,/^))0aM(^,/i) , (20) 

where X^a runs over all parton flavors. In Eq. (20), (j)a/A{^,f^) is the normal parton dis- 
tribution for finding a parton of flavor a in hadron A, and Ci/a = Z]n=o C*f"j(«s/vr)'^ are 
perturbatively calculable coefficient functions for finding a parton i in the modified parton 
distribution fi/A{xA,f^) from a parton a in the normal parton distribution 0a/^(^,/i). The 
first two coefficients of the Ci/a are available PJ26||: 



4°](^,^ = c/6) = 0, 



;i-^)+ --4U(1-;.) 



Cl/j{z,n = c/b) = 6ij— 

Cljl{z,fi = c/b)=TRz{l-z) , (21) 

where i and j represent quark or antiquark flavors and g represents a gluon. The coefficient 
functions given in Eq. (^) are derived from the general functional forms in Ref. |^ by setting 
the renormalization constants and the factorization scale as, Ci = c, C2 = 1, and fi = c/b. 

The (Tjj^y(M) in Eq. (0) is the lowest order cross section for a quark- antiquark pair to 
annihilate into the vector boson {V = 7*, W^, or Z^). For V = Z^, we have 

aij^zo{Mz) = Qij \ — ) , (22) 



where Qij is the weak charge. 



_ [l-4|eg|sin^g^y]^ + l 

with 9w denoting the Weinberg angle and denoting the charge of the quark. Here, i and 
j must represent a quark-antiquark pair of the same flavor. For the Weinberg angle we use 
sin^ Ow = 0.23. The aij^v{M) for V = W"^ and 7* can be found in Refs. @g. 



In the CSS resummation formalism, the Y term in Eq. (|I]) represents a small correction 
to the pt distribution when pt -C M. However, the Y term dominates the pr-distributions 
when Pt ~ M. The Y term has the perturbative expansion Y = Y,n=iY^^K'^sil-i') f'^)'^ , and 
the coefficients Y^"^^ have the factorized form [Bl 



t,A C,B 



where Y^a,b runs over all possible parton flavors and fi represents both the factorization and 
renormalization scales. The R^ab-^v ■'^l- (H) perturbatively calculable and have the 



same normalization as those introduced in Ref. |p. 

Since the perturbatively resummed Wij{h, M^xa^xb) in Eq. (|^) is only reliable for the 
small h region, an extrapolation to the large h region is necessary in order to complete the 
Fourier transform in Eq. @). In the original CSS formalism, a variable 6* and a nonpertur- 
bative function F^fgi^b, M,xa,xb) were introduced in such a way that W was modified 
for all b (except 6 = 0) as 

W'^'^ib, M, XA, xb) = W{k, M, XA, xb) Fc5s(^, M, xa, xb) , (25) 



where 6^, = b/ \ll + {b/bmaxY, with b^ax = 0.5 GeV^"*^. This construction ensures that 
6* < bmax for all values of b. The function F^gg is a Gaussian in 6, -F,^^ ~ exp(— /t6^) with 
K having a certain M^, x^, and xb dependence |0. 

As it has been pointed out in this paper and earlier |Ty], the original CSS extrapolation 
introduces a modification to the perturbative calculation, and the size of the modifications 
depends on the non-perturbative parameters in F^^g{b, M,xa,xb)- This motivated the 
development of the alternative extrapolation formula used in the present work. 
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FIGURES 
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FIG. 1. Integrand of the 6- integration in Eq. (|^ at = for as a function of h with an 
arbitrary normalization at = 350 GeV (dotted line), a/s = 500 GeV (dashed), and -y/s = 5.5 TeV 
(solid). 
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FIG. 3. (a) Cross section ^7 for production in pp collisions at RHIC with ^/s = 500 GeV; 
(b) defined in Eq. (|) for at the same energy. 
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FIG. 4. (a) Cross section for production in pp collisions at LHC with ^/s = 5.5 TeV; 
(b) defined in Eq. (S) for at the same energy. 
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FIG. 5. (a) Cross section for production in pp collisions at LHC with y/s = 5.5 TeV; 
(b) defined in Eq. (|) for at the same energy. 
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FIG. 6. (a) Cross section for production in pp collisions at LHC with y/s = 14 TeV 
defined in Eq. (^) for at the same energy. 
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FIG. 7. (a) Cross section for production in pp collisions at LHC with ^/s = 14 TeV; 
(b) defined in Eq. (|) for at the same energy. 
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FIG. 8. (a) The isospin ratio, Riso, defined in Eq. (|9|) (dashed hne), the shadowing ratio, Rsh, 
defined in Eq. (p!o|) (sohd hne), and "pure shadowing" (dotted) for Z*^ in p + Au cohisions at 
y/s = 350 GeV; (b) same as (a) for production. 
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FIG. 9. (a) The isospin ratio, Riso defined in Eq. (^) (dashed hne) and the shadowing ratio, 
Rsh, defined in Eq. (p^) (sohd line) for in Pb + Pb collisions at ^/s = 5.5 TeV; (b) same as (a) 
for production. 
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FIG. 10. (a) Cross section -^p^'^^''^'^\pT, Za/A, Zb/B) for production in p + Au collisions 
at RHIC with y/s = 350 GeV averaged over AB; (b) The ratio Rf^'^ defined in Eq. (|lT|) (soHd line) 
in p + Au collision and defined in Eq. (js]) for at the same energy. 
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FIG. 11. (a) Cross section ^'''^^''''\pt, Za/A, Zb/B) for production in Pb + Pb collisions 
at LHC with -y/s = 5.5 TeV, averaged over AB; (b) The ratio i?^^^ defined in Eq. (^) (sohd hne) 
in Pb + Pb collision and defined in Eq. (0) for at the same energy. 
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